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Abstract 

The strong-property-fluctuation theory (SPFT) provides a sophisticated means of estimating the 
effective constitutive parameters of a homogenized composite material (HCM), which takes account 
of the statistical distribution of the component particles. We present an extended version of the 
third-order SPFT in which the component particles are represented as depolarization regions of 
nonzero volume. Numerical results are provided for a bianisotropic homogenization scenario wherein 
the HCM is a Faraday chiral medium. Thereby, convergence of the extended SPFT at the second- 
order level of approximation is demonstrated within the long-wavelength regime. 
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1 Introduction 

Suppose that two (or more) homogeneous materials are blended together to form a composite. If 
the length scale of inhomogeneities in the composite is /3, then the composite is inhomogeneous 
at wavelengths < /3, but homogeneous at wavelengths ^ (5 (Lakhtakia, 1996). The homogenized 
composite material (HCM) which arises in the long-wavelength regime may exhibit properties that 
are not exhibited by its component materials, or at least not exhibited to the same extent. Thus, 
HCMs represent prime examples of metamaterials (Walser, 2003). Interest in complex metamate- 
rials, and HCM-based complex metamaterials in particular, has escalated in recent years, which 
serves to highlight the need for accurate formalisms to estimate the constitutive parameters of 
complex HCMs (Mackay, 2005). 

The strong-property-fluctuation theory (SPFT) provides a sophisticated basis for estimating the 
constitutive parameters of HCMs (Tsang &; Kong, 1981) which has distinct advantages over con- 
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ventional formalisms such as those named after Maxwell Gariictt and Bruggeman (Lakhtakia, 1996; 
Ward, 1995). This is achieved through accommodating higher-order descriptions of the distribu- 
tional statistics of the component materials which are brought together to form the HCM. Indeed, 
in principle, the SPFT can accommodate spatial correlation functions of arbitrarily high order. 
However, in practice, the SPFT is usually implemented at the second-order level of approximation 
wherein a two-point covariance function and its associated correlation length characterize the dis- 
tributional statistics of the component materials. Versions of the second-order SPFT have been 
developed for isotropic (Tsang & Kong, 1981; Michel &i Lakhtakia, 1995), anisotropic (Genchev, 
1992; Zhuck, 1994) and bianisotropic (Mackay, Lakhtakia, & Weiglhofer, 2000) linear HCMs, as 
well as for certain nonlinear HCMs (Lakhtakia, 2001; Mackay, Lakhtakia, &: Weiglhofer, 2003; 
Mackay, 2003). The third-order SPFT has also been established for bianisotropic HCMs which are 
weakly nonlinear. Convergence at the second-order level of approximation has been established 
for isotropic chiral mediums, and also more generally for bianisotropic mediums which are weakly 
anisotropic (Mackay, Lakhtakia, &; Weiglhofer, 2001a). 

Commonly, in homogenization formalisms — including the SPFT — the electromagnetic responses 
of the component material particles are represented by depolarization dyadics. Often the depolar- 
ization dyadics are taken to correspond to vanishingly small regions; accordingly, the spatial extent 
of the component material particles is neglected (Michel, 1997; Michel Sz Weiglhofer, 1997). An 
extended version of the second— order SPFT has recently been established which allows for de- 
polarization regions of nonzero volume (Cui Sz Mackay, 2007a). Numerical studies based on the 
extended SPFT have demonstrated that the depolarization contribution associated with regions of 
nonzero volume can have significant effects on estimates of the HCM constitutive parameters (Cui 
&; Mackay, 2007a, b) . We note that similarly extended versions of the Maxwell Garnett (Lakhtakia 
& Shanker, 1993; Shankcr &: Lakhtakia, 1993a, b) and Bruggeman (Prinkey, Lakhtakia, & Shanker, 
1994; Shanker, 1996) homogenization formalisms have also been established, and the spatial ex- 
tent of the component material particles has been emphasized in other homogenization studies too 
(Doyle, 1989; Dungey &; Bohren, 1991). However, these studies did not incorporate higher-order 
statistical details concerning the distribution of the component material particles, nor did they 
consider the most general linear scenario represented by bianisotropic HCMs. 

In this paper we present the extended third-order SPFT for biansotropic HCMs which are weakly 
anisotropic, and investigate the convergence of the extended SPFT. 

The following notation is used: Vector quantities are underlined. Double underlining and normal 
(bold) face signifies a 3x3 (6x6) dyadic. The inverse of a dyadic M is denoted by M~^. The 3x3 
(6 X 6) identity dyadic is represented by / ( I ) . All field-related quantities are implicitly functions of 
the angular frequency u. The permittivity and permeability of free space are denoted as e,, and /ig, 
respectively; the free-space wavenumber is ko = tu^eoHo . The real and imaginary parts of z G C are 
represented by Kez and Imz, respectively. A compact representation of the constitutive parameters 
for the homogeneous bianisotropic material specified by the Tellegen constitutive relations 
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is provided by the 6x6 constitutive dyadic 
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Herein, e and fi are the 3x3 permittivity and permeability dyadics , respectively, while ^ and 
( are the 3x3 magnetoelectric dyadics. Subscripts on K identify the particular material that K 
describes. 



2 Strong— property— fluctuation theory 
2.1 Component materials 

We consider the homogenization of two distinct material phases: phase a and phase b, both of 
which consist of spherical particles of average radius r/. All space is taken to be partitioned into 
the disjoint regions T4 and Vj, that contain the phases a and b, respectively. The phases a and b 
are randomly distributed, as specified by the characteristic functions 

r 1, reVi 

Mr) = { , ii = a,b). (3) 

I 0, r^Ve 

In particular, within the SPFT statistical moments of are utilized to characterize the component 
phase distributions. The volume fraction of phase i is given by the first moment; i.e., 

{^eir)) = fe, {i = a,by, (4) 

and we have fa + fb = ^- The physically-motivated step function (Tsang, Kong, &; Newton, 1982) 

r fe, |2:-r'|<L 

{Mr)Mr:)) = { , {e = a,b) (5) 

( fh \l-l'\> L 

is often adopted as the second moment for the second-order SPFT. The correlation length L is 
required to be much smaller than the electromagnetic wavelength(s) but larger than the size of 
the component phase particles. It is worth noting that the second-order SPFT estimates of the 
HCM constitutive parameters have been found to be largely insensitive to the particular form of the 
second moment (Mackay, Lakhtakia, &: Weiglhofer, 2001b). In keeping with (5), the third-order 
SPFT has been established for the third moment (Mackay, Lakhtakia, & Weiglhofer, 2001a) 

ff, min{Li2,Li3,L23} > L 



{Ml) Ml') Ml")) 



fe, max{Li2, -Li3, L23} < L 

{£ = a,b), 

l{fi + 2ff), oneoiLi2,Li3,L23<L 
I {2fe + ff), two of L12, Lis, L23 < L 

(6) 



3 



where 



(7) 



The component material phases a and h are homogeneous materials, characterized by the 6x6 
constitutive dyadics and K^, respectively. 



2.2 Homogenized composite material 



The constitutive dyadic of the HCM, as estimated by the nth-order SPFT, is given by (Mackay, 
Lakhtakia, & Weiglhofer, 2000) 

1 r . , 1-1 



kI"] = K 

= HCM =cm 



I + St"J(?7,L) -0(7?) 



S^"l(r/,L). 



(8) 



Herein, the constitutive dyadic characterizes a comparison medium whose constitutive pa- 

rameters are provided by the Bruggeman homogenization formalism. 

The depolarization dyadic 

Uiv) =U° + U^\v) (9) 

has two parts (Cui Sz Mackay, 2007a): D° represents the contribution to the depolarization arising 
from component particles in the limit ry — ^ 0, whereas 0^*^(7/) represents the depolarization con- 
tribution arising from the nonzero volume of the component particles. Often in homogenization 
studies the 0^^(17) contribution is neglected, but recent studies have highlighted the significance 
of this contribution, particularly in the context of scattering losses (Cui & Mackay, 2007a,b). The 
conventional SPFT incorporates D° only as the depolarization dyadic, whereas the extended SPFT 
accommodates both D'' and 0^^(77). The mathematical expressions for D° and 0^*^(17) are com- 
plicated, especially for bianisotropic HCMs, but integral representations are available which can 
be straightforwardly evaluated using standard numerical techniques (Press, Flannery, Teukolsky, 
&L Vetterling, 1992). These integral representations are provided in the Appendix. 

The mass operator term S^"' (ry, L) in (8) vanishes for the zeroth- and first-order versions of the 
SPFT (Tsang & Kong, 19"81; Mackay, Lakhtakia, & Weiglhofer, 2000); i.e., 

^[0] ^ ^[1] ^ g_ (^Q) 

By implementing the two-point covariance function (5), the second-order mass operator term is 
given by (Mackay, Lakhtakia, & Weiglhofer, 2000) 



SPl(r?,L) 



fafb X (.rj) -X(v) 
= a =0 



D>0(L) 



X iv) -x(v) 

= a =0 



with the polarizability density dyadics 



(11) 



{e = a,b). (12) 



The three-point covariance function (6) yields the third-order mass operator term (Mackay, Lakhtakia, 
k Weiglhofer, 2001a) 



HI -fay =a' 



Y{V)'X (r/)-D>°(L) 



+ D>0(L) 'X iv)' Yiv) + D>"(L) • X iv)- D>"(L) • x iv), (13) 
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where 



Y(r?) = -K;^-D(r/). 



(14) 



3 Numerical studies 



We now apply the extended third-order SPFT presented in §2 to a specific bianisotropic ho- 
mogenization scenario. As an illustrative example, let us consider the homogenization of (i) a 
magnetically-biased ferrite medium described by the constitutive dyadic (Lax & Button, 1962; 
Collin, 1966) 
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and (ii) an isotropic chiral medium described by the constitutive dyadic (Lakhtakia, 1994) 
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The parameter 5 in (15) provides a means of varying the constitutive contrast between the com- 
ponent material phases. The constitutive relations of the resulting HCM — which is known as a 
Faraday chiral medium — are rigorously established (Engheta, Jaggard, &; Kowarz, 1992; Weigl- 
hofer & Lakhtakia, 1998). The constitutive dyadic of the HCM, as estimated by the nth order 
SPFT, has the general form 
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(17) 

An HCM of the same form also arises from the homogenization of a magnetically-biased plasma and 
an isotropic chiral medium (Weiglhofer, Lakhtakia, & Michel, 1998; Weiglhofer &; Mackay, 2000). 
Bearing in mind that the third order SPFT is established only for bianisotropic mediums which are 
weakly anisotropic (Mackay, Lakhtakia, & Weiglhofer, 2001a), we select the following representative 
values for the constitutive parameters of the component material phases: = 1.2 + i0.02, /j,^ = 
3.5 + ii0.08, fii = O.T+ iO.OOS, fi^ = 3.0 + 7:0.06; = 2.5 + iO.l, = l + i0.07 and fib = 1.75 + ?;0.09. 
Results which are qualitatively similar to those presented here were observed — in further studies 
not reported here — when different values were selected for the constitutive parameters of the 
component materials. 
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In the following numerical studies, the correlation length L is fixed^ for each value of 5 considered, 
while the particle size parameter r] varies from to L/2. In order to conform to the long-wavelength 
regime under which the SPFT estimates of the HCM parameters are derived, the value of L is 
selected such that the scalar Q <^1, where 

Q = ma^{l7i|, I72I, I73U74I} ^^g^ 
27r 

with I i = 1, . . . , 4} being the four independent wavenumbers supported by the HCM. For sim- 
plicity, we choose the wavenumbers associated with propagation along the Cartesian z axis (Mackay 
& Lakhtakia, 2004); i.e.. 



71 — ^0 yy^HCM + ^HCMy ^^^HCM + ^''HCM ^HCM ^HCM 



C9 



(19) 



73 = ^0 y\l ^%CM ~ '^HCM\j ^liCM ~ i^HCM + ^HCM ~ ^HCM 

74 = K (y-\J^^HCM ~ ^%CM\J f^HCM ~ f'HCM + ^HCM ~ ^HCm) 

All numerical calculations were carried out using an angular frequency a; = 27r x 10^° rad s~^ with 
the volume fraction fixed at /„ = 0.3. 

The zeroth-, second- and third-order SPFT estimates of the HCM constitutive parameters I^hcm^ 
^^^HCM ^'^'^ l-''HCM ^^'^ plotted against the size parameter r] in Figure 1 for the case where 5 = 10. 
Here, the correlation length is set at L = 0.45 mm, in order that Q = 0.1 at 77 = L/2. The 
third-order estimates of the real and imaginary parts of /J-h^m, I^hcm I^hcm increase steadily 
as 77 is increased, as do the zeroth-order and second-order estimates, but the difference between 
the third-order estimates and the second-order estimates remains very small for all values of ij. 
In contrast, there are plainly significant differences between the second-order and zeroth-order 
estimates. Furthermore, the difference between the zeroth- and second-order estimates increases 
in magnitude slightly as the the size parameter 77 increases. The corresponding graphs for the 
permittivity and magnetoelectric constitutive parameters of the HCM are qualitatively similar to 
those graphs presented in Figure 1. Accordingly, these are not displayed here. 

Plots of iJ-^Qjij, f^HCM ^^hcM versus 77 for the case where 5 = 30 are shown in Figure 2. The 
correlation length L = 0.28 mm was used for the calculations of Figure 2, thereby resulting in 
Q = 0.1 at 77 = L/2. As is the case for 5 = 10, the second-order estimates are plainly different 
to the zeroth-order estimates for 5 = 30. The differences between second-order and third-order 
estimates of the real and imaginary parts of the HCM constitutive parameters increase slightly 
as 5 increases, but they remain small for all values of 77. The corresponding permittivity and 
magnetoelectric constitutive parameters of the HCM exhibit trends which are qualitatively similar 
to those shown in Figure 2 for the HCM magnetic constitutive parameters. 



^This applies to the second-order and third-order SPFT calculations; L does not feature in the zeroth-order 
SPFT. 
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4 Concluding remarks 



It is demonstrated by our numerical studies in §3 (and in further studies not presented here) that 
the extended SPFT at the third-order level of approximation does not add significantly to the 
HCM estimates yielded by the second-order extended SPFT. The differences between second- 
order and third-order estimates of the HCM constitutive parameters arc very small for all values 
of the size parameter investigated, even when the constitutive contrast between the component 
materials is as large as a factor of 30. Significant differences between the second- and third- 
order estimates arise only when (i) the correlation length and/or size parameter become similar 
in magnitude to the electromagnetic wavelength(s); and/or (ii) the constitutive contrast between 
the component materials becomes enormous. In the case of (i) the bounds imposed by the long- 
wavelcngth regime are exceeded, while in the case of (ii) the contrast between the polarizability 
density dyadics % and is likely to become strong. In either scenario the basic assumptions 

underlying the long-wavelength SPFT are violated (Mackay, Lakhtakia, & Weiglhofer, 2000, 2001a). 
We therefore conclude that the extended SPFT converges at the second-order level of approximation 
for bianisotropic HCMs which are weakly anisotropic. 
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Fellowship. 



Appendix 

The depolarization dyadic D(r/), as specified by (9), for a spherical inclusion of radius rj, immersed 
in a bianisotropic comparison medium described by the constitutive dyadic K^^^, is derived from 
the dyadic Green function of the comparison medium. The contribution associated with 77 — >■ is 
provided by the ?7-independent surface integral (Michel, 1997; Michel &; Weiglhofer, 1997) 




(20) 



with dyadic integrand 



B {q) = 



1 r Oi^{6,(l))qq 



-ac_{e,<j))qq 



(21) 



b{^A) [ -a^{9,(p)qq 



Herein the scalars 



ap{e,(j))=q'P 



= cm 



(P = e,C,C,/x) 



(22) 



and 



b{e, 4>) = [a,{e, <p) a^ie, <^)] - [a^ie, cp) a^ie, </>)] , 



(23) 
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while q is the radial unit vector specified by the spherical coordinates 9 and (f>. 

The depolarization contribution associated with the nonzero inclusion volume is given by the rj- 
dependent volume integral (Mackay, 2004; Cui & Mackay, 2007a) 



D>0(r?) 



27r2 



ILU 



sm(qr]) 
QV 



— cos{qr]) 



(q) - B (q) 



d^q, 



(24) 



wherein 



A (g) 

= cm ^— ' 

als (20) a 

techniques (Press, Flannery, Teukolsky, & Vetterling, 1992). 



+ K 



(25) 



g {q/u) X I 

-{q/uj) X X g 

The depolarization integrals (20) and (24) are straightforwardly evaluated using standard numerical 
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Figure 1: Real (left) and imaginary (right) parts of the HCM constitutive parameters fJ'^cM plotted 
against rj (mm) for 6 = 10. Key: dashed curve is the zeroth-order SPFT estimate; broken dashed 
curve is the second-order SPFT estimate; and solid curve is the third-order SPFT estimate. 
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Figure 2: As Figure 1 but for 6 = 30. 
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